Kondo Effects in Quantum Dots at Large Bias 
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Recently, the issue of whether the Kondo problem in quantum dots at large bias is a weak-coupling 
problem or not has been raised. In this paper, we revisit this problem by carefully analyzing a 
corresponding model in the solvable limit — the Emery-Kivelson line, where various crossover energy 
scales can be easily identified. We then try to extract the scaling behavior of this problem from 
various physical correlation functions within the spirit of "poor man's scaling." Our conclusions 
support some recent suggestions made by Coleman et al. [Phys. Rev. Lett. 86, 4088 (2001)], 
which are obtained by perturbative analysis: The voltage acts as a cutoff of the renormalization 
group flow for only half of the impurity so that the low-temperature physics is controlled by a strong- 
coupling fixed point. But the low-temperature response functions in general show one-channel Kondo 
behaviors with two-channel Kondo behaviors occurring only through proximity to a quantum critical 
point. 

PACS numbers: 73.63.Kv, 72.10.Fk, 72.15.Qm 



I. INTRODUCTION 

The Kondo problem is one of the best studied many- 
body problems in condensed matter physics. Due to ad- 
vances of nanotechnology, the Kondo effect in quantum 
dot systems predicted by theoriescl was observed in recent 
experiments.13 Although the quantum dot is intrinsically 
a multilevel system, as the energy is much lower than the 
single-particle level spacing, the system can be described 
by a single-impurity Anderson model with the level spac- 
ing playing the role of cutoff in this model. Moreover, 
in the Coulomb blockade regime with an odd number of 
electrons, it can be mapped ontO|-a two-lead Kondo model 
with the following Hamiltonian:tl 
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momentum k and spin a, and Jl, Jr, and Jlr(= Jrl) 
are positive (antiferromagnetic) Kondo coupling con- 
stants between electrons and the dot (S). Energies 
e Q k = £k — eV a , where V a = ±V/2 are the potentials of 
the left- and right-hand leads. Besides, derived from ail 
Anderson model via a Schrieffer- Wolff transformations, 
the coupling constants obey \Jlr\ 2 = JlJr- 

Within the spirit of "floor man's scaling," a recent pa- 
per by Kaminski et al.p shows that even in the present 
nonequilibrium case, the low-temperature properties can 
still be characterized by a single crossover energy scale 
Tk which is identified as the Kondo temperature of this 
problem. Through extending the information gained 
from perturbative renormalization group (RG) equations, 



Kaminski et al., use the ordinary one-channel Kondo 
fixed point to extract low-temperature transport prop- 
erties for V < T K . However, for V > T K , the RG- 
improved perturbative calculation indicates that the con- 
ductance will saturate at a value much smaller than 
that in the unitary limit. This signals that the cou- 
pling between different leads Jlr stops growing for en- 
ergy smaller than V . One then wonders whether the 
above-mentioned one-channel Kondo (1CK) fixed point 
can still be used to describe the low-temperature physics 
in the latter case. Recently, to-,gain further insight into 
this problem, Coleman et al.p have done a perturba- 
tive calculation of the impurity magnetic susceptibility. 
The scaling behaviors for various couplings they obtained 
are basically the same as that obtained in Ref. || in the 
high-temperature region. However, for T <C V , the flows 
°f Jl.r and Jlr exhibit different behaviors: Jlr stops 
growing at the energy scale V while Jl,r continue to 
grow toward strong coupling with a new Kondo temper- 
ature T K . They conjecture that the physics at T < T* K 
is described by a 2CK fixed point. (See also Ref. ||.) 

In this paper, we shall revisit this problem by studying 
the model in s. solvable limit — the generalized Emery- 
Kivelson line.H This enables us to follow the RG flow 
all the way from the perturbative regime down to the 
low-temperature region, and the concept of universal- 
ity guarantees that we can obtain correct scaling at low 
temperature. Moreover, through studying various cor- 
relation functions in details, we can see how the impu- 
rity is screened (or unscreened) due to the presence of 
finite bias. This approach is complementary to the per- 
turbative analysis in Ref. 4, which is supposed to be 
valid at high temperature. Combined with the results 
obtained from the perturbation theory, we reach a com- 
plete picture about the behavior of the Hamiltonian (jl]) 
at large bias: The effect of the cotunneling term H tra ns is 
to generate a new energy scale Tlr even in the channel- 
symmetric case {Jl — Jr)- As Tlr <§C Tk, there exists 
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a range of temperature T^r <C T <C 2V, where the uni- 
form magnetic susceptibility exhibits the 2CK behavior, 
i.e., the logarithmic temperature dependence. On the 
contrary, it will show the 1CK behavior as Flr rs Tr. 
This is the fundamental difference between the present 
case and the ordinary 2CK,fixed point. Previous studies 
on ordinary 2CK problemsQ revealed that the presence of 
an unscreened Majorana fermion lies at the heart of the 
2CK properties, which is the origin of the logarithmic 
temperature dependence in various suscepetibilities. In 
our case, the logarithmic temperature dependence is due 
to the partial screening of a Majorana fermion. Although 
the low-temperature behavior of the uniform susceptibil- 
ity is similar to that in the ordinary channel-asymmetric 
2CK problem ( Jl ^ Jr. and H tr ans — 0) and there are 
two characteristic energy scales in both cases, the origins 
are distinct. In the ordinary channel-asymmetric 2CK 
probelm, the new energy scale in addition to the Kondo 
temperature is generated from the channel-asymmetry, 
i.e., Jl Jr, whereas in our case it arises from Htrans 
and still exists even J L = J R . 

The rest of the paper is organized as follows. We in- 
troduce the solvable model in Sec. II. In Sec. Ill, the im- 
purity Green functions and the impurity contributions to 
the uniform magnetic susceptibility are calculated. The 
last section is devoted to a discussion and conclusions of 
our results. 



II. THE SOLVABLE MODEL 

After changing the notation rS(ik) — * 1(2), we start 
with the following Hamiltonian:t!lli3 
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Hq + Hi + H2 



where 
H 
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H 2 = 



H.c. 



(2) 



Here a = 1,2, a =T,|, and V 2 = -Vi = V/2. For 
simplicity, we set e = 1 and the Fermi velocity Vp = 1 
and : • • • : is the normal ordering with respect to the 
Fermi surface in the absence of V. (Here we assume that 
Fermi velocities on both leads are equal.) 

The procedure to arrive at a solvable model is given 
in the following. First we have to bosonize the Hamilto- 
nian (0) . To deal with Klein factors carefully, we employ 



bosonization formulas on a finite length 
1 



s 



F aa exp{-iA L (N a<7 - P Q /2)x} 



x exp {— iV4^4>aa} 



(3) 



where N aa is the number of ip aa fermions and Al = 
2n/L. Here L is the system size and ao is a short- 
distance cutoff. Pq takes care of the boundary conditions 
of fermion fields. Fa^s are Klein factors, which satisfy 
the commutation relations [F aa , N a '<r'] — 5aa'5 aa >F aa , 
= 2<W<W, and {E aa , F a , a ,} = 0. Next 
four bosonic fields are introduced] by 
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and the corresponding transformation on N aa is 
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(5) 



Here N m e Z + P/2 for m = c, s, /, sf and £ m = 
mod 2. P = 0, 1 for the total number of fermions 
being even and odd integers, respectively. After plugging 
these into the Hamiltonian (Q), we perform the Emery- 
Kivelson (EK) transformation: U = exp {i^/4irS z (t) s (0)}. 
To proceed, we also introduce four more Klein factors 
F m with m = c,s,f,sf, which satisfy the commutation 
relation: [F m ,N m r] = S mm 'F m . With the help of Eq. 
(||), the identification between F aa and F m can be found. 
What we need is the following ones: F^Fi^ = F s F s f, 

F^F 2 i = F.F« f , F 2 \F 1T = FfF s , and F^F n = F f F*. 
Then, the EK-transformed Hamiltonian is refermionized 
by the following formulas: 



$ = F S S+, d=FjS-, 
S z = dU-l/2, 



(6) 



and 
V m (x) 

* s (x) 



1 



F m exp {-iA L (N m - l/2)x} 



V27TO0 

x exp{-iV47T0 m }, m = /, sf, 

' F s e l7rdfd exp{-iA L {S T - l/2)x} 



x exp{-i\/47r0 s } ! 



(7) 
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with St = S z + N s being the total spin operator. Finally, 
the structure of the Hamiltonian can be further simplified 
by introducing Majorana fermions 



V2 
d + dt 
y/2 



y/2i 



, rn= /, sf, 



d-d) 
s/2i 



(8) 



Now, by taking L — > oo and ignoring terms of 0(\/L) , 
we obtain 

H' = UHU ] =H a + H f + H sf + H int + const, 
where 

dx : ^lid x ^ s (x) 



Hn 



1 b \c b ic 



q>0 



H 



iV^Ifl^(O), 



Hi, 



ab {i6X 1 :^ s (0):-8X 2 e sf e sf (0) 

+ 8X 3 e sf e f (0)} ■ (9) 



Here b nr — 



qc 

Af x /(47rao) with : 



y/2q/L J dxe iqx <j> c {x) (q > 0). = 
i = ±, LR. SXi = A+|| — 2-7T , <5A2 = A_y , 
and <5A 3 = A Lfl ||. X ±l = (X u ± A 2i )/2 with i =_L, ||. No- 
tice that for symmetric dots, T_ = = SX 2 - The charge 
sector is completely decoupled from the impurities and 
we will not consider it hereafter. From Eq. (|9|), we see 
that the EK line corresponds to SXi = (i = 1, 2, 3) and 
its Hamiltonian Hf + H s f can be solved exactly. 



III. PHYSICAL OBSERVABLES 

Now we are in a position to compute impurity con- 
tributions to the uniform magnetic susceptibility, which 
is defined by \%mp = 

\\m L ^oo(dM / dh\ h=ti - Lxo)/L 
where \o is the Pauli susceptibility of bulk electrons 
and M = <7(Sr) is the magnetization. (We consider the 
case with gyromagnetic ratios on the impurity site and 
bulk electrons being equal and denote it by g.) From 
Eq. ([?]) and the operator product expansion (OPE) 
of *|(z + ia)^ s (z), we have S T = J dx : *|.* s (a;) : 
+ (P + l)/2. The last term does not depend on the mag- 
netic field h and we ignore it. 

In the presence of external magnetic fields, there is 
an additional term in the Hamiltonian, AH — —ghSx, 
which is invariant against the EK transformation. Per- 
forming the EK transformation and refermionization suc- 
cessively, we obtain the Hamiltonian describing noninter- 
acting fy s fermions as H s = J dx : ^\{id x — gh) 1 i' s (x) : . 



The term linear in h can be removed by the transforma- 
tion $ 8 (x) -► # s (x)e~ l9hx . As a result, H s -> f dx : 
^lid x ^ s (x) : +const.. By noticing that : *|* s (a;) :->: 
\&J\& s (a;) : +gh/(2ir), the magnetization becomes 



M = g / dx(: *t* s (x) :) + L Xo h, 



(10) 



while Hi n t turns into Hi nt = Hi nt + i(gh/2ir)5Xiab. It is 
straightforward to see that on the EK line the first term 
in Eq. ([h]) vanishes. This leads to Ximp = and the 
leading contribution to Ximp must arise from Hi nt .u 

To calculate contributions to Ximp given by Hi nt , we 
use the Keldysh formula.EJ Before plunging into the cal- 
culation, we need impurity Green functions G a (ui) and 
Gb(ui) on the EK line: 



G a (u) 

G b (u) 











-2i 



r I , J? .(n v _ a; -n v + ^)+r_(l-2n u) ) 

uj 2 +r 2 



-iT 



(11) 



u—ir^ 



where T = T^r + T_ and n u is the Fermi distribution 
function. Here all Green functions are defined in the 
Keldysh space as 



Gr Gk 
G A 



G = 



An important feature in Eq. ( Jl 1| ) is that in the presence 
of the cotunnelling term {Tlr 7^ 0) a fermions are always 
screened by electrons in leads at low temperature (<C T) 
even for symmetric dots. However, we shall see later 
that the interpretation off depends on the ratio |V|/T+. 
Finally, the bare Green functions of ^ s fermions are 



G>, P ) = 







2n t 
i 



,)5{uj+p) 



u)-\-p— i0+ 



(12) 



With the help of Eqs. fly) and (|2J) and H int , we 
are able to find that in the perturbative expansion of 
magnetization, the leading nonvanishing order is the one 
with SXi- It consists of two terms: one independent of h 
while the other proportional to h. The former becomes 
zero as V — and it does not contribute to the magnetic 
susceptibility. Consequently, the latter gives the leading 
behavior of Ximp at low temperature. After performing 
integrals over p and u> and in terms of digamma function 
ip{x), we have 



Xirnp(T, V) 



(gsx.f 



i 



4?r 3 r + - r 



with 



h 

h 



+ 



1 

2 ' 2ttT 
T LR 



-if, 



^Re {it - , 

r+ + r r V 2 2ttT 



iV 



2 + 2ttT J ' 

r + + iV 



1 2 2ttT 



(13) 
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Here Re{/(z)} means the real part of f(z). I\ gives the 
result of the ordinary 2CK problem and T + plays the role 
of Kondo temperature Tr- in that case without channel 
asymmetry. Effects of finite bias on Ximp are through 
the function I2 arising from scattering in the spin-flip 
cotunneling channel. 

Based on Eq. and in terms of the asymp- 

totic formula of diagamma function ip(z) = 
lnz — l/2z — l/12z 2 + • • •, we can discuss the low 
temperature behavior of Ximp- It depends on the ratio 

R = r+/T. 

(i) R = 0(1). When T < T + ,f , we have 

h = in(r + /f) + o([T/r + (f)] 2 ), 

h ~ \nf(V)+0{(T/e 1 f), (14) 



with f(V) = ^(T\ + V 2 )/(T 2 + V 2 ). Here the crossover 

energy e x = T+(f) and \V\ for \V\ < T+(f) and \V\ > 
r + (r), respectively. Thus, the leading behavior of Ximp 
is 



Xrmp(T, V) = 



(9SX1) 2 



1 



4tt 3 r+ - r 



4 



■lnf(V) 



(ii) R > 1. When r < T < T+, we have 



47T 3 



In (T£/T), 



(15) 



(16) 



with Tpr 



t k {i + [r L fl/2(r+ + r)] in [i + (r+/v) 2 ]} 

and T K for \V\ > T and \V\ < T, respectively. Here 
A = 1, 1 + r Lfl /(r + + f ) for \V\ > T and |V| < T, 
respectively. Tr- = cr + with c = 2d 1 /it and 7 being the 
Euler constant. Since Tlr/T + <C 1, Tr- is in general 
with the same order as Tk- For extremely low tempera- 
ture T <C r, the behavior of Ximp turns into that shown 
by Eq. dl5|). It is, however, a crossover not a phase tran- 
sition. Equations (15) and ( |l6| ) are the main results of 
this paper. We shall discuss their implications now. 



IV. CONCLUSIONS AND DISCUSSIONS 



Combined with perturbative analysis in Ref. |I| we ar- 
rive at the following picture: In the case with r + « T, 
all couplings flow to the strong-coupling regime at low 
energy as |V| <C T + and both a and b fermions are com- 
pletely screened. It leads to 1CK behavior. This can also 
be understood by noticing that the two-lead Hamiltonian 
Hf + H s f in Eq. (||) has the structure of two copies of 
the 2CK problem where half of the impurity is screened 
by ipf and another half of the impurity is screened by 
■0s/, and the voltage acts only on ipt. (For simplic- 
ity, we take T„ = 0.) For V —> and T < T LR ,T + , 



both channels will flow to strong couplings and the sys- 
tem manages itself into a single-channel Kondo prob- 
lem, which can be easily seen by taking T + = Tl,r 
and rewriting the Hamiltonian through defining a new 
fermion tp — (1/V2)(£? + *£]/)■ On the other hand, for 
\V\ ^S> T + , the flow of Jlr stops at the scale V while 
Jl,r still continue to flow towards strong couplings with- 
out being affected by the voltage. This is the origin of 
the scale e\ appearing in Eq. (jlj). It also leads to the 
fact that the Majorana fermion a cannot be completely 
screened and fluctuates with a scale f . Therefore, we 
expect that the 2CK behavior emerges in some situation 
as shown in Eq. (|l6|). This in part explains the stability 
of this 2CK problem against the perturbation Jlr- It 
is, however, not exactly equivalent to an ordinary 2CK 
problem because the coupling Jlr below the scale V is 
not irrelevant but marginal. In other words, half of the 
impurity, the a fermion, can be viewed as a free fermion 
only at energy scales higher than T. As the temperature 
is far below it, the impurity behaves as if it is completely 
screened. 

The partially screened Majoiana fermions also reveal 
themselves in the conductance^ 



(17) 



where 4>'{z) — dip/dz and Gjj — (Tlr/T 1 )e 2 /ir is the 
conductance in the unitary limit. The low-temperature 
asymptotic behavior for T <C T is 



G(T, V) 



G, 



Y 2 + V 2 



0((T/e 2 ) 



(18) 



with the crossover energy e 2 = \V\ and f for | V\ ^> f and 
|V| <C f , respectively. From Eq. (|l8|), we see that the 
conductance exhibits very different behaviors for large 
and small bias. For bias much smaller than the crossover 
energy T, the conductance will reach the unitary limit, 
while for V ^> T, the voltage plays the role of a crossover 
energy scale, and the conductance will saturate at a much 
smaller value than that in the unitary limit. Within the 
spirit of "poor man's scaling," we can extract scaling be- 
haviors of the corresponding coupling: For V <tiT, Jlr 
will flow towards strong coupling and completely screen 
the impurity. However, for 7 > f, the RG flow will 
terminate at the scale V and the impurity can not be 
completely screened. Note that the scaling behavior of 
Jlr gained here is completely consistent with that ob- 
tained in Ref. |]. In addition, the asymptotic behavior 
of Eq. ([l8]) is similar to that appearing in Ref. ||, which 
reveals the existence of some kind of universality in this 
nonequilibrium problem. 

Concerning the possible experimental realization of 
the above-mentioned 2CK behavior, we need to be close 
enough to that quantum critical point. This requires a 
large ratio R for the renormalized scales, which is equiv- 
alent to the condition Jlr <C Jl.r- However, it is very 
difficult to achieve this goal due to the constraint on the 
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bare couplings, J\ R ~ JlJr, and the logarithmic nature 
of the RG flow. Thus, unless we can design a dot with a 
very small ratio for bare interlead and intralead tunnel- 
ings, it seems that this window is too narrow to observing 
the 2CK behavior. 

Finally, we would like to mention that similar crossover 
behaviors, i.e., the 2CK behavior at the intermediate- 
temperature regime while the ordinary 1CK behavior 
in the extremely low-temperature regime, have already 
appeared jp the context of channel-asymmetric 2CK 
problems.E3 Although some technical details are not com- 
pletely the same, the underlying physical reasons are sim- 
ilar — the appearance of a new energy scale suppresses 
some scattering processes and changes the direction of 
the RG flow. We would, also like to point out that a re- 
cent paper by Zvyaginliil discussed a similar problem but 
in a totally different physical context by using the Bethe 
ansatz. More precisely, Ref. |ll| considers the low-energy 
properties of conduction electrons hybridized with local- 
ized 5/ electrons. When the concentration of 5/ elec- 
trons is low, the magnetic susceptibility exhibits a sim- 
ilar crossover behavior to our Eqs. ( |l5| ) and (|l6|). In 
that case, the new energy scale in addition to the Kondo 
temperature arises from the hybridization anisotropy. In 
a word, the most crucial difference between the models 
considered in Refs. [T(] and |fl| and the present one is that 
the appearance of two crossover scales in the former case 
is due to channel anisotropy, while in the latter case it 
is due to the cotunneling processes. The mechanisms to 
generate the new energy scale reveal their distinctions in 



the channel-symmetric limit. In that limit, for the mod- 
els considered in Refs. flOl and 11, the new crossover scale 



vanishes and the low-temperature dynamics is described 
by the 2CK fixed point. On the contrary, for the Hamilto- 
nian ([!]) in the channel-symmetric limit, the two crossover 
energy scales never vanish, and depending on the relative 
magnitude of the two scales (a large difference induced by 
large bias), the low-temperature dynamics is controlled 
either by the 2CK or 1CK fixed points. However, it al- 
ways exhibits lCK-type behaviors when the temperature 
is far below both scales. 

To sum up, we arrive at the following conclusions: 
(i) Regardless of the magnitude of bias V, the low- 
temperature physics of quantum dots coupled to two 
leads is controlled by a strong- coupling fixed point, (ii) 
This strong-coupling fixed point should exhibit the be- 
havior of a one- channel Kondo fixed point instead of a 
two-channel one. (iii) There is a nearby two-channel 
Kondo fixed point located at the unphysical parameter 
space f = ( in the sense that T = requires Jlr = 0). 
For dots with T <C Tk, it controls the physics at the range 
of temperature I <T« Tk ■ 
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